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ABSTRACT 

Pulsations driven by partial ionization of hydrogen in the envelope are often considered important for driving 
winds from red supergiants (RSGs). In particular, it has been suggested by some authors that the pulsation 
growth rate in a RSG can be high enough to trigger an unusually strong wind (or a super-wind), when the 
luminosity to mass ratio becomes sufficiently large. Using both hydrostatic and hydrodynamic stellar evolution 
models with initial masses ranging from 15 to 40 M0, we investigate 1) how the pulsation growth rate depends 
on the global parameters of supergiant stars, and 2) what would be the consequences of a pulsation-driven 
super-wind, if it occurred, for the late stages of massive star evolution. We suggest that such a super-wind 
history would be marked by a runaway increase, followed by a sudden decrease, of the winds mass loss rate. 
The impact on the late evolution of massive stars would be substantial, with stars losing a huge fraction of 
their H-envelope even with a significantly lower initial mass than previously predicted. This might explain the 
observed lack of Type II-P supernova progenitors having initial mass higher than about 1 7 M0 . We also discuss 
possible implications for a subset of Type Iln supernovae. 

Subject headings: stars: evolution, stars: mass-loss, stars: massive, supergiants, supernovae:general 



1. INTRODUCTION 

Mass loss due to stellar winds is one of the governing fac- 
tors that determine the evolution of massive stars. Great im- 
provements in our understanding of line-driven winds from 
massive hot stars (i.e., O-Type stars or Wolf-Rayet stars) 
have been achieved over the last three decades, and they 
are well reflecte d in recent stellar evolution models (see 
iPuls et alTI I2OO8I for a recent review). However, the driv- 
ing mechanism for cool-star winds remains elusive. Note 
that many of the authors stUl re ly on the empirical mass- 
loss rate given by Ide Jager et al. I (1988) (hereafter, JNH88), 
which is more than 20 years old, for modeling cool giant stars, 
while some alternati ve prescriptions have been considered 
by so r ne authors (e.g.lSalasnich et al. II 19991: Ivan Loon et al. I 
120051: IVanbeverenetal. 1 120071) . Late stages of massive star 
evolution during and after core helium burning are therefore 
relatively difficult to study, which is a major obstacle to theo- 
retical investigation on supernova progenitors and their envi- 
ronments. 

One of the key factors that should be taken into account 
for the study of the red supergiant (RSG) evolution is pulsa- 
tion. It is well known that RSGs are unstable to radial pul- 
sations driven by partial ionizati on of hydrogen in the outer- 
most layers of the envelope (e.g. lLi & Gong II 19941) . and pul- 
sation is one of the most commo nly invo ked m echanisms for 
driving winds from cool stars ( cf. lBowen 1988; H ofner et al. I 
120031: iNeilson & Lester1l2008l) . Although the detailed mech- 
anisms for driving mass loss by pulsations are not well un- 
derstood, simple arguments based on the gain of pulsation 
energy suggest that RSGs might experience a stronger mass 
loss for a high er growth rate of the pulsational instabili ty (e.g. 
lAppenzeUaf|[T970.) . In particular. fHeger et al. I ( 119971) found 
that RSG pulsations can achieve large amplitudes that may be 
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comparable to those in AGB stars when the pulsation peri- 
ods become close to the thermal time scale of the envelope. 
They therefore suggested that such strong pulsations might 
lead to a "superwind" if the luminosity to mass ratio in a RSG 
were sufficiently high. This idea might b e supported by obser- 
vations. For example, Ivan Loon et al. I {2005) found that the 
dust-enshrouded RSGs in the Large Magellanic Cloud have a 
very high mass loss rate, which appear s to be related to strong 
RSG pulsation (Ivan Loon et al. II2008I) . 

In the present study, we explore possible consequences of 
such a superwind in the evolution of supernova progenitors. 
Our approach to this topic is as follows. First, we investi- 
gate the properties of RSG pulsations using our new stellar 
evolution models that follow the non-linear evolution of pul- 
sation (Sect. |2]i. This allows us to derive a relation between 
the pulsation growth rate and the structure of RSGs. Second, 
we assume that the mass loss rate due to RSG winds should 
be significantly enhanced compared to that of JNH88, when 
the growth rate of pulsation becomes sufficiently large. The 
consequent evolution of massive stars can be very different 
from the standard picture given in previous work (Sect. [3]). 
Note that here we apply the mod ified mass loss prescrip tion 
only to very unstable RSG s, while lSalasnich et akl (Il999h and 
IVanbeveren et al. I (l2007h considered overall enhancement of 
RSG mass loss rate for all types of RSGs. Also, our mass loss 
prescription is more directly related to the pulsational prop- 
erties of RSGs than the others that rely on surface luminosity 
and temperature. Finally, we conclude the paper by discussing 
implications for supernova progenitors (Sect.|4]i. 

2. PROPERTIES OF RSG PULSATIONS 

RSG pulsations are not usually found in stellar evolution 
models even if one uses a hydrodynamic code, for two rea- 
sons. First, in most stellar evolution calculations, the evolu- 
tionary time steps are much larger than the pulsation periods 
in RSGs, which are comparable to the dynamical time scale 
of the envelope. Second, pulsations tend to be numerically 
damped out with the implicit method adopted in stellar evolu- 
tion codes. However, the study bv.He ger et al. , (.1997.) shows 
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Fig. 1. — Evolutionary tracks of hydrostatic reference models with initial 
masses of 17, 20, 25, 30 and 40 M0, on the HR diagram. The color shad- 
ing gives pulsation periods (upper panel) and growth rates (lower panel; see 
the text for the definition), obtained with the hydrodynamic calculations that 
are performed at the given reference points. The open diamonds mark the 
reference points where no pulsation is detected with our code. 



that the predictions of linear stability analysis for basic pul- 
sational properties (i.e., pulsation period and growth rate) can 
be accurately reproduced in non-linear evolutionary calcula- 
tions of RSGs by a hydrodynamic stellar evolution code, if the 
growth rate is sufficiently large and if the adopted time steps 
are sufficiently small. On the other hand, non-linear evolu- 
tionary calculations do not require the condition of complete 
thermal and hydrostatic equilibrium of the star, in contrast to 
the case of linear stability analysis. We therefore decide to 
use an existing hydrodynamic stellar evolution code for our 
analysis of RSG pulsations . Details of our im plicit hydrody- 
namic code are described in lYoon et al. I(^2006^ and references 
therein. 

To start with, we construct non-rotating hydrostatic models 
at solar metallicity for 11 different initial masses (15, 16, 17, 
18, 19, 20, 21, 23, 25, 30, 40 Mq). Since the hydrodynamic 
term in the momentum equation is switched off, pulsations do 
not appear in these model sequences. The HR diagrams for 
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Fig. 2. — Radial velocity at the suii'ace in a hydrodynamic RSG model 
sequence with Mi„it = 20 Mq . The surface luminosity and the effective tem- 
perature at the reference point are logL/hQ = 5.05 and T^ff = 3198 K, respec- 
tively. 
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Fig. 3. — Pulsation growth rate versus some physical quantities of RSGs. 
The solid line gives the best fit, as expressed in Eq. (T). See the text for more 
details. 

some model sequences are shown in Fig.[T] 

Then, we turn on the hydrodynamic term and restart the 
calculation from the reference models chosen at various evo- 
lutionary epochs during the supergiant phase for each model 
sequence. In total, 71 reference models are selected as start- 
ing points. In Fig. [1] some of the starting points are shown, 
marked by filled circles and open diamonds. To guarantee 
enough time resolution for describing pulsations in the hy- 
drodynamic calculations, we limit the time steps such that 
At < O.Ol[/?/(lOOOR0)]2[(15Mo)/M] yr. This choice is 
based on the fact that pulsations of a RSG with Mjnjt > 15 Mp, 
have a period of the order of 1000 days (e.g. iHeger et al. I 
119971) . and obey the relation of P oc R^/M, where P is the 
pulsation period, R t he radius and M the total mass of a RSG 
(IGoughetal.1[T96l . 

An example of such calculations is given in Fig. |2] As 
the pulsation amplitude grows gradually, the surface veloc- 
ity eventually reaches the sound speed. Since our code cannot 
describe shock waves, we investigate the pulsation properties 
only with the sub-sonic results. Fig.[T]shows the pulsation pe- 
riod P and the growth rate ry of the fundamental mode at dif- 
ferent reference points on the HR diagram. Here, the growth 
rate is defined as 77 = \v{to + P)/v(to)\, where v(fo) denotes a 
local maximum of the surface velocity located at a certain 
time t = to- According to this definition, pulsational instability 
means r/ > 1 . 
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The growth rate obtained here should be considered only 
indicative, given the limitations of our adopted numerical 
method. For example, in our calculations, no pulsation is de- 
tected from some reference points. This does not mean that 
the considered RSGs are stable against the pulsational insta- 
bility, but that the growth rate is too small for our code to 
follow. Also, in our calculations, instant adjustment of con- 
vective flux during pulsation is assumed, which is not phys- 
ical given the similar time scales of convection and pulsa- 
tiorU Despite these uncertainties, our results are qualitatively 
in good agreement with those of previous studies: the growth 
rate is generally higher for a larger luminosity to mass ratio 
(L/M), and/or for a smaller thermal time scale of the enve- 
lope. Specifically, as shown in Fig. [3] we find the follow- 
ing relation between the pulsation growth rate and the stellar 
structure from our results with the hydrodynamic model se- 
quences: 



log?7 = Ci 



R 



M 



'nCH,er 

yr 



-0.315 



-Cz, (1) 



where TKH.env (:= GMiotMenv / LR) denotes the Kelvin- 
Helmoltz time scale of the hydrogen-rich envelope, Ci = 
9.219 X 10"'' and C2 = 0.0392844. 

We confirm that RSGs are generally more susceptible to the 
pulsation instability for larger initial masses. Fig. [T]indicates 
that for Minit > 20 M0, pulsations can become strong enough 
to be detected with our code even before the star reaches the 
Hayashi line, while pulsations are only found after neon burn- 
ing for Minit = 17 Mq. We do not find any signature of pul- 
sation for Minit =15 M(r|3- Therefore, if a super-wind could 
be induced by strong pulsations, it would occur in an earlier 
evolutionary epoch for a higher mass star. 

3. THE EVOLUTION OF SUPERGIANTS WITH SUPER- WINDS 

If the mass loss rate were dramatically enhanced due to 
strong pulsations, how would it affect the evolution of super- 
giants, and what would be the observational consequences? 

To address this issue in more detail, we perform an exper- 
iment as the following. We assume that the mass loss rate of 
a supergaint is significantly enhanced compared to the rate of 
JNH88, if pulsations are strong enough to be detected with our 
hydrodynamic stellar evolution code, i.e. if 77 > 1 in Eq. ([T]i: 



M=r/"MjNH88, 



(2) 



where Mjnhss is the mass loss rate of JNH88, and a a free pa- 
rameter. New evolutionary model sequences are constructed, 
with this modified mass loss rate. Here, the hydrodynamic 
term is switched off again, and normal evolutionary time steps 
are used. 

It should be kept in mind that the empirical mass loss rate of 
RSGs of JNH88 might already reflect the role of pulsations. 
However, as implied by the growth rates shown in Fig.lT] only 
very unstable massive stars with Minit ^19 M© are signif- 
icantly affected by this modification. Given the paucity of 
RSGs w ith Minit >20 - 25 M^ in our galaxy and nearby galax- 
ies (e.g. iMassev et al.1l2009h . our assumption that the mass 



^^ Note, however, non-linear calculations bv'He ger et al. 1 1 19971) give quite 
consistent results with those of their linear stability analysis that assumes a 
frozen-in convective flux 

"* Heger et al. ( 1997) found pulsations for the same initial mass. The rea- 
son for this difference is that their models consider rotationally enhanced 
mass loss, resulting in a higher L/M than in our non-rotating models. 
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Fig. 4. — Upper panel: The evolution of 20 and 25 M© stellar models 
on the HR diagram. The thick solid lines and the thin dashed lines give 
the evolutionary tracks with a = 1.75 (i.e., with pulsation-induced mass loss 
enhancement), and a = 0.0 (i.e., no enhancement of mass loss) in Eq. (2). 
respectively. The calculation was terminated at a point close to the end of 
helium burning (central helium mass fraction equal to 0.0013) for the 25 Mq 
model sequence with a = 1.75, while the evolution is followed up to neon 
burning for the other sequences. The filled circle on each track with a = 1.75 
marks the position when the mass loss rate reaches its peak. The filled star 
gives the position of the pre-supernova star for the 20 Mq sequence a = 1.75. 
Lower panel: The mass loss rates of the corresponding stellar models. 

loss history of such very unstable RSGs can be quahtatively 
different from that of most RSGs in the sample of JNH88 does 
not necessarily cause an inconsistency. 

Fig.|4]shows the evolutionary tracks and the mass loss rates 
of model sequences with a= 1 .75 (this choice for a is rather 
arbitrary), compared to the corresponding ones with a = 0, for 
two different masses (Minit = 20 and 25 Mq). Two important 
points are worthy of notice here. 

First, as the star gradually loses mass, the L/M ratio and the 
pulsation growth rate increase accordingly. If the value of a 
is large enough (as adopted here), this leads to a runaway in- 
crease in the mass loss rate as shown in the lower panel of the 
figure. Second, the mass loss rate decreases very rapidly when 
the star moves away from the pulsationally unstable regime 
on the HR diagram. The peak of the mass loss rate is reached 
when the helium mass fraction at the centre decreases to 0.54 
(0.031), and the total mass to 9.35 M© (6.50 M©), for 25 Mq 
(20 M0) model sequence. As shown for the 20 M© sequence, 
such a mass loss history of runaway increase followed by sud- 
den decrease can be repeated as the star moves in and out the 
unstable regime. The 20 M© star finally ends its life as a yel- 
low supergiant of Mtot = 6.1 M©, with only a small amount 
of hydrogen of about 0.5 M©. The most likely outcome of 
the death of such a star would be a Type lib supernova. The 
25 M© star has a hydrogen envelope of 0.22 M© at core he- 
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Hum exhaustion, and would produce either a Type lb or a lib 
supernova depending on the subsequent history of mass loss. 
When the mass loss rate reaches the maximum, the growth 
rate given by Eq. ([U becomes as high as 1 1.3 and 10.6 for 25 
and 20 Mq models, respectively. We checked if non-linear 
evolutionary calculations also give such high growth rates, 
given that the relation of Eq. Q is only based on the result 
with t; < 8.0 (Fig. |3]l. We find that, at the reference points 
marked by the filled circles in Fig.Hl the hydrodynamic calcu- 
lations give 77 = 9.9 and 9.0 for 25 and 20 Mq models respec- 
tively, which are comparable to the values given by Eq. ([U. 

4. DISCUSSION 

If a pulsationally-driven super-wind (PDSW) phase could 
be induced by strong pulsations, this would have very impor- 
tant implications for supernova progenitors. As implied by 
our model sequences presented above, such a PDSW can sig- 
nificantly reduce the critical ZAMS mass (Mdit) above which 
a huge fraction of the stellar envelope is removed before core 
collapse, thus producing no Type II-P supernova. Based on 
our results, we suggest Md-it ~ 19...2OM0 with a PDSW (see 
Figs.[T]and|4]i, while the models with the JNH88 mass loss rate 
predict SN II-P progenitors up to ~ 25M0. Interestingly, this 
value of Merit is comparable to what othe r alternative prescrip- 
tions of RSG winds mas s loss predict ( Salasnich et al. I[l999t 
IVanbeveren et al. 120071) . Note that Merit could b e even further 
reduced by rotation, as implied by the result of lHegeret al. I 
( II99 7) (see also Meynet & Maeder 2003). 

This might provide a plausible solution to the so-called 
RSG problem, i.e. the observe d lack of type II-P progen- 
itors with Mi„it > 16.5 ± 1 ■5MQ (ISmartt et al. I l2009h . Such 
observation could also result from the presence of dusty cir- 
cumstellar material around the stars, as an obscured progen- 
itor would be estima ted to have a lower initial mass in pre- 
SN images. However ISmartt et al. I (.200 9) also noted that the 
number of SNe type II-P in their sample is consistent with the 
expected number of stars in the range 8.5-16.5 Mq assum- 
ing a Salpeter IMF. This could be a possible indication of the 
fact that stars with higher initial mass do actually lose most of 
their hydrogen envelope, which cannot be easily understood 
with the canonical mass loss rate of JNH88. 

The PDSW is expected to significantly affect the circum- 
stellar medium around a RSG (van Veelen, in prep.). When 
the star dies, the shock produced by the collision between the 
SN ejecta and the circumstellar material transforms kinetic 
energy into thermal energy. This energy can be radiated away 
at different wavelengths, res ulting in brightn ening the super- 
nova remnant for long times (IChevalier 1 1 977h . If the collision 
occurs directly after the SN explosion, the emission can even 
alter the spectrum and light curve of the SN. This scenario 
is usually invoked to explain the sub-class of Type Iln super- 
novae (Schlegel 1990; Filippenko 1997) . To reproduce the 
light curve of the most luminous Type Iln supernovae (SNe 



Iln) like SN 2006gy and SN 2006tf, very massive shells are 
required, which need to be eje cted in eruptive events a few 
years before core collapse (e.g., 'van Marie et al. II2OIOI) . Pul- 
sational pair instability (Wooslev et al. 2007f) and LBV-like 
eruptions (e.g., 1 Smith et al. 1 12007!) have been discussed to 
explain the presence of shells of 1O-25M0 around the pro- 
genitors of luminous SNe Iln. A PDSW phase is an unlikely 
explanation for such extreme environments. However it might 
be interesting for those type Iln where the required mass in the 
stellar vicinity is of the order of a few solar masses or less. 

In this context it is interesting to consider the circumstellar 
medium around the RSG VY CMa. The stellar surrounding 
appears shaped by episo dic mass ejections which occurred 
about 500-1000 yr ago ( Smith etal.1 120091) . The complex 
morphology of the CSM is suggestive of a possible interaction 
between convection and pulsation (see Fig. 13 in Smith et alT] 
2009), which for these stars ar e predicted to occur on similar 
timescales (e.g., iHeger et al. I 1997). The mass loss rate of 
1 -2 X 10"^ Meyr"' derived by Smith et al. (2009) is compa- 
rable to the one expected during the PDSW phase (see Fig.|4li. 
Even if the mass loss prescription for PDSW we used is some- 
what arbitrary, we want to stress that the energy available from 
the growth of pulsa tion is enough to dr ive mass loss rates up 
to '^ 10"^ M0yr"'. ISmith et al~l (l2009h state that an exti-eme 
RSG like VY CMa would produce a Type Iln event like SN 
1988Z if it were to explode in its current state. Therefore the 
occurrence of pulsation-driven super- winds in RSG might ex- 
plain moderately luminous SN Iln if the enhanced mass loss 
takes place less than about ^ 10^ years before core collapse. 
With the mass loss prescription of Eq. |2] with a sufficiently 
large a (^ 2), this would occur in a narrow range around 
^19M0 for non-rotating stars. An accurate determination of 
the expected rate of type Iln due to PDSW requires a more 
realistic mass loss prescription. 

To conclude, if a pulsationally-driven super-wind phase 
could be induced by strong pulsations, we would expect a sub- 
stantial change in the late evolution of single massive stars. 
The mass loss rate during the RSG phase would increase dra- 
matically for stars with Mi„it > Merit compared to the JNH88 
rate, and the PDSW phase should start earlier in the evo- 
lution of more massive stars. The resulting pre-supernova 
structure of these stars is affected, as well as their CSM. For 
single stars, this suggests the following sequence in super- 
nova types as function of increasing initial mass: II-Ph> II- 
L— ^ Iln— )► lib— )> lb— > Ic. Detailed estimate for the mass range 
of each supernova requires a better prescription for the pulsa- 
tionally enhanced mass loss, which deserves further study. 

We are grateful to Stan Woosley, Bob van Veelen and Nor- 
bert Langer for many useful discussions. S.C.Y is supported 
by the DOE SciDAC Program (DOE DE-FC02-06ER41438). 
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